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Abstract 

A D-like meson resonance with mass around 2.9 GeV has been found in the DKK system 
using two independent and different model calculations based on: (1) QCD sum rules and (2) 
solution of Fadeev equations with input interactions obtained from effective field theories built by 
considering both chiral and heavy quark symmetries. The QCD sum rules have been used to study 
(— | . the -D s *o(2317)-^ and D/o(980) molecular currents. A resonance of mass 2.877 GeV is found with 

the D/o(980) current. Although a state in the D s * o(2317)-fr current is also obtained, with mass 



around 2.9 GeV, the coupling of this state is found to be two times weaker than the one formed in 
D/o(980). On the other hand, few-body equations are solved for the DKK system and its coupled 
channels with the input t-matrices obtained by solving Bethe-Salpeter equations for the DK, DK 
and KK subsystems. In this study a D-like meson with mass 2.890 GeV and full width ~ 55 
MeV is found to get dynamically generated when DKK gets reorganized as -D/o(980). However, 
no clear signal appears for the D s *q(2317)K configuration. The striking similarity between the 
results obtained in the two different models indicates strongly towards the existence of a D/o(980) 
molecule with mass nearly 2.9 GeV. 
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I. INTRODUCTION 



In the past years, the development of high energy facilities has lead to the discovery 
of a number of open and hidden charm resonances by collaborations like BABAR, Belle 
and BES 1 |j] which, in turn, has motivated many theoretical studies to understand the 
properties and nature of heavy flavor hadrons. Among which some of the heavily discussed 
states are D s * (2317), X(3872), Z + (4430), whose properties have been studied within dif- 
ferent models assuming different configurations like diquarks, tetraquarks, hybrids, hadron 
molecules, etc. (for a review see Refs. [6|-l9|). 

The understanding of the nature of the different mesons and baryons of the hadron 
spectra, in general, is a long standing puzzle in theoretical nuclear physics. QCD is the 
accepted fundamental theory describing the strong interactions in terms of the quarks and 
gluons which constitute the hadronic matter. However, while at high energies the theory 
becomes perturbative and has been successfully tested by the experiment, the situation is 
very different at low energies, where due to the confinement of the quarks the theory is 
not anymore perturbative and nonperturbative methods are needed to extract information 
about the properties of the hadrons. 

To face this challenging issue different techniques have been developed. One of them 
is Lattice QCD, which in the last few years has emerged as an important tool to extract 
information about hadronic observables like mass, phase shifts, etc. However, due to the 
large number of degrees of freedom present in QCD (quarks and gluons of different flavors 
and colors) numerical calculations involving big number of lattice points and small lattice 
spacing are very time consuming for natural values of the mass of the quarks. Although a 
lot of progress has been done in this area, there are still some problems when addressing 
excited states which have decay channels lot 14 1. 

Another a 
sum rules 
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temative to study hadrons within the spirit of QCD is the method of QCD 
. In this formalism the hadrons are described in terms of their interpo- 



lating quark currents, with which a correlation function is built. One begins evaluating 
this correlation function at short distances, where the quark-gluon dynamics is essentially 
perturbative, and then nonperturbative corrections are added to it. This method has been 
widely used to understand the mass, coupling, decay width, etc., of many hadron states. 
Yet another way to elucidate the nature and properties of mesons and baryons is based on 



the use of effective field theories built by taking into account unitarity, chiral symmetry and 
its spontaneous breaking. In this case, the hadrons are the degrees of freedom of the theory 
instead of the quarks which constitute them. In the last 20 years, there has been lot of ac- 
tivity in this field and many resonances have been found to have important meson-meson or 
meson-baryon components in their wave functions. Some of the states most widely discussed 
are the A (1405), generated as a consequence of the interaction of the coupled channel system 



KN and vrS 
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j|-|23j . and the /o(980) resonance, formed in the KK and irir system 24j-|26|. 
Recently, this theory has been generalized to study the properties of hadronic systems in 
a finite volume and its value in the determination of related physical observables using the 
energy levels obtained in the finite volume and, thus, as a prospective tool for Lattice QCD 



calculations, has been shown 



23-l21|. 



The above mentioned methods are in continuos development since the experimental ac- 
cess to higher and higher energies is becoming plausible and, consequently, more and more 
new states with heavy quarks are being found. Present time is thus ideal to study heavy 
hadron physics since model predictions can be immediately tested, which eventually helps 
in understanding the structure of hadrons. With this idea we present a study of the DKK 
system, which we find particularly interesting since the DK and KK interactions are at- 
tractive in nature. In this manuscript, we have studied this system using two methods: 
QCD sum rules and Few-body equations. In the former case, we investigate the Dj o(980) 
and D s *o(2317)K configurations, while in the latter we solve the Faddeev equations for the 
three-hadron system, where /o(980) and -D<j*o(2317) are dynamically generated in the cor- 
responding subsystems. As we shall see, we find a resonance with similar characteristics in 
both models. 

In the following, we first discuss the calculations based on QCD sum rules and the results 
found in it. Subsequently, we tackle with the formalism to solve the Faddeev equations and 
discuss the results obtained with it. Finally we draw some conclusions. 
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II. QCD SUM RULES 



We start our study based on the QCD sum rules by writing the interpolating molecular 
currents for the Df and D S * K systems as 

jOfo = i fe75Ca ) ^ b s b ) (1) 

f**o R = i (s a c a ) (q bl5 q b ) , (2) 

where a and b are color indices, and q represents a light quark (u or d). Using these currents, 
we write the two-point correlation function 

I% 2 ) = z j d*xe^(0 | T [j(x)j\0)} | 0), (3) 

which can be written in terms of the quark propagators by contracting all the quark anti- 
quark pairs (for more details see, for example, Ref [6|). 

This function is of a dual nature: it represents a quark- ant iquark fluctuation for short 
distances (or large negative q 2 ) and can be treated in perturbative QCD, while at large 
distances it can be related to hadronic observables. The sum rule calculations are based 
on the assumption that in some range of q 2 both descriptions are equivalent. One, thus, 
proceeds by calculating Eq. ([3]) for both cases and by eventually equating them to obtain 
information on the properties of the hadrons. 

From the QCD side, for large momentum transfers, Eq. fl3]) can be calculated, in the 
first approximation, by assuming the involved propagators as those of free quarks. However, 
since we are finally interested in studying the properties of hadrons, the relevant energies 
are lower, where the distance between the quarks gets longer and quark-gluon interactions, 
quark-antiquark pair creation becomes important. It is, thus, required to include the effect of 
the presence of the gluons and quarks in the QCD vacuum. For practical calculations, then, 
one resorts to the Wilson operator product expansion (OPE) method, where the correlation 
function is expanded in a series of local operators 

n OPE = ^c n (Q 2 )o n . (4) 

n 

In Eq. 01]) the set {O n } contains all local gauge invariant operators expressible in terms of 
the gluon fields and the fields of light quarks and the coefficients C n (Q 2 )(Q 2 = —q 2 ), by 
construction, include only the short-distance domain and can, therefore, be evaluated per- 
turbatively. Nonperturbative long-distance effects are contained only in the local operators. 
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In the expansion of Eq. (J4]), the operators are ordered according to their dimension n, 
where n = corresponds to the unit operator, i.e., perturbative contribution, and the 
rest of operators are related to the QCD vacuum fields in terms of condensates. For normal 
quark-antiquark states, the contributions of condensates with dimension higher than four are 
suppressed by large powers of Aq CD /Q 2 , with 1/Aq C d being the typical long-distance scale. 
However, for molecular states, condensates with higher dimensions can play an important 
role. This is taken into account by writing Eq. (jl]) in terms of the spectral density using the 
dispersion relation 

oo 

r n OPE (s) 

n OPE / 2\ = / dg h + Subtraction terms, (5) 

\ j s — q 2 

ml 

with 

p° PE (g 2 ) = p pert + p ma + p^ ^ + p ms ^ -\- p(i9<*-Gq) _|_ pin s (qgcr-Gq) _|_p<OT> 2 _|_ p™* <99> 2 ? f§\ 

and m c representing the mass of the charm quark. The spectral density p OPE is related to 
the imaginary part of the correlation function as 7rp OPE (s) = Im [n OPE (s)]. 

A calculation of the different terms in Eq. ([6]) for the D s *qK and Dfo systems leads to 
the following expressions, where m s is the mass of the strange quark. 

1. The perturbative or dimension contribution is found to be: 

.pert / 2% «ert, 2% _ f , ((<* ~ l)g +mj a 



2. For the terms of dimension 1, which are proportional to m s , we get: 

r _ {q2) = _ da m c ((a-l)q +m c ) a m s 



D s *qK ' J 2 10 7T 6 (a- l) 3 



C (g 2 ) = 0. (8) 

3. The calculation of the diagrams with one quark condensate gives 

3m c ((a — l)g 2 + m 2 ) 2 a 2 (s s) 

2 ' 



{qq) f q 2\ = / d 

H D s , K Kq ' / 2% 4 (a - 1 





OLmax 



P D f0 [q) - J ^ 2%4( a _l)2 ' W 
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4. Both p ms (ii) and p^ G% ) contribute to dimension four and the expressions for the 
corresponding spectral densities are, 

,2 i ™2\ 2 



n mMi) {( ?\- f , 3((a-l)g + m c ) am s 
p D 3 , R [q >- J m 2^\a-l) 



(qq) - (ss) 



nia x 



p m a{qq){ 2 ) = _ I 9 ((a - l)q 2 + ml) 2 a(s s}m s 

P D fo [q ' J 2W(a-l) ' 1 J 

o 



p(^ 2 _) (g 2 )=p(9 ^> (g 2 )= / da 3 (9 2 G 2 ) / (2-a)(m c 2 + (q-l)g 2 ) | m 2 c a 2 



d„*oR ' ' Df Q ^ 1 J 2 12 tt 6 V 2(1 - a) 9(1 - a) 2 







w a(mg + {a - l)g 2 ) 

X 1^ (Uj 

5. Considering the mixed condensates, we get 

3m c (sga ■ Gs)a(l — 2a){m 2 c + (a — l)g 2 ) 



p { ™ a - Gq) {q 2 ) = / da 



r>,. * ^ y J 2 8 vr 4 (a - 1) 
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Jqg*-G q ) fn 2\ f A „ 3m c (qga ■ Gq)a(l - 2a)(m 2 c + (a - l)q 2 ) 

P D f0 [q) ~ J da 2H\a-l) 2 [ ] 

o 

6. Going to the dimension 6 operator, we get the following contributions for the terms 
proportional to m s {qga ■ Gq), 

OLrnax 



CT G * >( " 2)= / <*"sfiW-« , < 1 -» 





(qga ■ Gq) (6Ln (a) — 3) — (sga ■ Gs) ( — — 

V 1 — a 



m s (sga ■ Gs) 



P™f^V)= J da '"'^ w/ [ml - q 2 (1 - a)] (1 - 6Ln (a)) , (13) 
o 

two quark condensates 

{qq? ( 2 ) = _ da ^~^ +m)(qq)(ss) 

D s * aK 1 J 2% 2 



^V)= /■fa ((a - 1) ^ m ' )< " ) , (14) 
o 
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and three-gluon condensates 



7. Finally, for dimension 7, we get 

c| V)=/ da ^(^-( ? - ? )(M 

o 

p-<*«>V) = - / rfa^Mf^ (16) 

'iJ/o W ^ J 2% 2 V 7 
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The integration limit in Eqs.(l71)-(|T6l) is a max = 1 — ^f-. For numerical calculations we need 
the values of the different condensates and quark masses, which are given in Table [D 

TABLE I: Values of the different known parameters required for numerical calculations of 



the correlation function given by Eq. fl3]) (see Ref. 32|) 



Parameters Values 

m s 0.10 ± 0.022 GeV 

m c 1.23 ± 0.05 GeV 

{qq) -(0.23 ± 0.03) 3 GeV 3 

(ss) 0.8 (qq) 

(g 2 G 2 ) 0.88 GeV 4 

(g 3 G 3 ) 0.045 GeV 6 

(qa ■ Gq) 0.8(gg) GeV 2 

We now calculate the correlation function from the hadronic or phenomenological point 
of view. In this case, the currents j* and j are interpreted as the creation and annihilation 
operator of the hadrons which have the quantum numbers of the current j. For this n(g 2 ) 
is written by inserting a complete set of states with the same quantum numbers as those of 
the currents under consideration 

n pheno m(g2) = • f d 4 xeiq . X f d 3 P ! j{x) ! m ^ (m ^| jt(()) | 0) _ 

J J 2p°(27r) ^ 



Thus, the correlation function contains the information on all the hadrons of a given set of 
quantum numbers including the one we are interested in, which is the low mass, relatively 
narrow, hadron of the series. One proceeds in such a situation by assuming that the spectral 
density of hadrons, for a fixed set of quantum numbers, can be expressed as a sum of a 
narrow, sharp state (which we are interested in), and a smooth continuum 

p phcnom (s) = X 2 6{s _ ^ + ( s ) , (1 8 ) 

where p C ontmuum is assumed to vanish below a certain value of s, so, which corresponds 
to the continuum threshold. Above this threshold, it is assumed to be given by the result 
obtained with the OPE. Therefore, one uses the ansatz |33j Pcontinuum(s) = p OPE (s)0(s — s ). 



In Eq. ( 1181) A is the coupling of the current j with the low-lying hadron with mass m: 
(0\j\m) = X. 

The spectral density given by Eq. (fTST) is related to the correlation function of Eq. f[T7|) 

as 

oo 

n pheno m(g2) = + / ds P W (lg) 

m — Q J s — Q 

so 

To carry out the calculations, So is taken as a parameter of the method but its value is 
not completely arbitrary: it is related to the onset of the continuum in the current j under 
consideration and is taken to be roughly 0.5 GeV above the mass of the hadron we are 
interested in jf], [l^. In this work, we are looking for a resonance with a possible Df o(980) 
or D s *o(2317)K molecule-like structure. Since such resonances are weakly bound, they are 
expected to get generated close to the threshold of the constituent mesons. Thus, a/s in 
the present case can be ~ 3.4 GeV. 

The correlation function calculated using QCD suffers from divergent contributions com- 
ing from long range interactions, while the one calculated phenomenologically contains con- 
tribution from the continuum. This situation can be improved by taking Borel transform 
of both Eqs. flH]) and (fi~9l) . which kills the problematic terms of both sides. After taking 
the Borel transform, we equate the resulting expressions of the correlation functions on the 
basis of its dual nature and get 

oo oo 

A 2 e -mV^ + J dsp OPE (s)e- s / M2 = j dsp OPE (s)e- s / M \ (20) 

so rr& 
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which can be rearranged as 

SO 

x 2 e -mVAP = j dsp OPE( s)e -s/M^ (21) 

ml 

where M represent the Borel mass parameter. Calculating the derivative of Eq. (1211) with 
respect to M 2 and dividing the resulting expression by Eq. (I2ip . we obtain the mass sum 
rule 

J ds s p OPE (s)e- s / M2 

m* = =L . (22) 

/ dsp OPE {s)e- s / M2 

ml 

Having the mass one can evaluate the current-state coupling constant through Eq. (l2~Tj) 

/ dssp OPE (s)e~ s / M2 
^ = ~ e _ m 2 /A /2 • (23) 

The reliability of the results obtained within QCD sum rules depends on the definition 
of a valid Borel window. This range of the Borel mass is obtained by making the following 
constraints: 

• The maximum value of the Borel mass, M max , where the results should be reliable, is 
fixed by ensuring that the pole term (low mass hadron) gives the dominant contribution 
to the calculations. However, M max is a function of sq. As mentioned earlier, a 
reasonable value of y/s in the present calculation can be 3.4 GeV. We show the 
contributions of the pole and continuum terms weighted by their sum {(], [y, 35] for 
the Df {980) and D S * (2317)K systems, obtained with ^ = 3.4 GeV, in Fig. HJ 
which shows that M^ ax ~ 2.05 GeV 2 in the former case and 1.77 GeV 2 in the latter 
one, respectively. 

• The second constraint is to look for that Borel mass range where a convergence in the 
OPE series is found. For this we calculate the perturbative contribution and add to 
it the diagrams with higher dimensions step by step. In other words, we calculate the 
right hand side of Eq. (12 lj) by first using Eq. (ED) for p OPE , then by using the sum of 
Eqs. Q and ([8]), which means including the diagrams up to dimension 1, next we do 
the calculations up to the subsequent higher dimension by taking a sum of Eqs. ©-(19]), 

9 




c ■ 

o r . i ... i ... i ... i ... i ... i ... i ... i ... i . i r . i ... i ... i ... i ... i ... i ... i ... i ... i . 

° 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 

M 2 (GeV 2 ) M 2 (GeV 2 ) 



FIG. 1: The contributions of the pole (solid line) and continuum (dashed line) weighted by 
(divided by) their sum for the D/ o (980) (left panel) and D S * (2317)K (right panel) 

systems. 



etc, until going to diagrams with dimension 7 (given by Eqs. (Tl6l) ). For a convenient 
comparison, the result obtained in each case is weighted (divided) by the one obtained 
by using the whole series of Eq. fl6]) for the spectral density. In Fig. |2j we show the 
results of such an analysis of OPE convergence for the D/o(980) (left panel) as well 
as D S * (2317)K (right panel) systems. 



Df (980) 




1.8 

M 2 (GeV 2x 



D S . (2317)K 



^=3.4 GeV 



Dim 
Dim 1 
Dim 3 
Dim 4 
-Dim 5 



I I — — Dim 6 



1.5 



1.6 



1.7 

M 2 (GeV 2 ) 



1.9 



FIG. 2: Relative contributions of the different OPE terms as a function of the squared 
Borel Mass, for the -D/o(980) (left panel) as well as D s *q(2317)K (right panel) systems. A 
value of \/s = 3.4 GeV is used in these calculations. The arrows in the figures indicate the 

valid Borel window, which is determined by using the conditions discussed in the text. 
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The final condition, which is imposed to identify the minimum value for the Borel 
mass, is that the contribution defined by 



E [fdsp°™(s) 



s/M 2 



E [fdsp°^(s)e 

dim=l \ m 2 



s/M 2 



Nmax / so 



dim=l 



(24) 



is less than 0.25. In the equation written above, Nmax refers to the maximum di- 
mension of the condensates taken into account in the calculation, which is 7 in the 
present case. With this, the OPE convergence is ensured in the Borel mass range 
where the results can be taken as the reliable ones. It can be see from Fig. [2] that a 
good convergence of the OPE series is found for M Min = 1.5 GeV 2 and 1.74 GeV 2 , for 
the -D/o(980) (left panel) and _D s *o(2317).K" (right panel) cases respectively 

The valid Borel windows established using both criteria discussed above are indicated 
with arrows for the ,0/0(980) and D S * (2317)K systems in Fig. El for ^ = 3.4 GeV. 

Having fixed these conditions, we show, in Fig. |3l the results of the calculation of the 
mass (using Eq. ( 1221) ) of the states described using Df (980) and D S * (2317)K molecule-like 
currents. It can be seen that a mass of 2.867 ± 0.004 GeV is found in the case of _D s *o(2317)-K" 



D f (980) 




D S . (2317)K 



M 2 (GeV 2 ) 




M 2 (GeV 2N 



FIG. 3: Mass of a molecule-like resonance obtained by using QCD sum-rules to study the 
-D/o(980) (left panel) and D S * Q (2317)K (right panel) systems. The arrows indicate the 
range where the reliability of the results is ensured. 



while it is 2.925 ± 0.032 GeV in the Df (980) system (within the Borel window, indicated 
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by the arrows in Fig. [3]). However, these calculations have been done using the value of 3.4 
GeV for y/s , which is a parameter. We now vary the continuum threshold in the range 
3.3 < y/sQ < 3.6 GeV to check the sensitivity of our results to this parameter. This leads to 
a mass for the resonance in D s *o(2317)K system as 2.963 ± 0.1 GeV, while in the case of 
£>/o(980) we get 2.98 ± 0.1 GeV. 

We have also checked the sensitivity of the mass to the uncertainty in the knowledge of 
the values of the different condensates and quark masses listed in Table. HI Considering all 
these uncertainties we finally get 

m Ds , oK = (2.978 ± 0.15) GeV, (25) 

for the D S * (2317)K molecular current, while for the Df (980) molecular current we get 

m Dfo = (2.877 ± 0.36) GeV. (26) 

Next, we have calculated the coupling A for the two configurations studied and found 
that the current-state coupling of the state described by the D s *o(2317)^ current is around 
two times weaker than the one found for the Df (980) current: 

\ DstoK = (6.4 ± 1.8) x 10~ 3 GeV 5 , and A D/o = (9.3 ± 5.3) x 10~ 3 GeV 5 . (27) 

We can interpret this result as an indication that a Df (980) state is better represented by 
the respective molecular current than the D s *o(23l7)K state. 



III. THREE-HADRON APPROACH 



Let us now discuss the study of the DKK system within a very different approach which 
is based on effective field theories, treating hadrons as the degrees of freedom instead of 
quarks, and examine if the findings obtained in such a calculation are compatible with the 
ones found with QCD sum rules. 

l36l . l37| . the /o(980) resonance is described as a 



In the unitary chiral models [19 



-26 



molecular hadron state generated in the interaction of the KK, and tttt coupled channels 24 , 
251 ] . Similarly, the -D s *o(2317) state can be interpreted as a DK bound state formed in 



the DK, D s t] coupled channel system 3814421] . In these models, Lagrangians based on 



symmetries like chiral 



43N45| and heavy quark symmetries 46j-|48j are used to determine 
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the lowest order amplitude describing the transition between the different coupled channels. 
These amplitudes are further unitarized by using them as driving terms in the Bethe-Salpeter 
equation, and the scattering matrix t for the system is obtained. Recently, these models 
based on effective field theories, chiral symmetry and unitarity in coupled channels have 
been further extended to investigate the interaction of three-hadron systems formed by 
different mesons and baryons, like, irKN, NKK, J/ipKK, <j)KK, etc., and generation of 
several hadron states like £(1660), 0(2170), F(4260), iV*(1710), has been found |4jH32]. 
Analogously to the two-body studies where the scattering matrix is obtained by solving the 
Bethe-Salpeter equation taking as kernel the lowest-order chiral amplitude, in the case of 
the approach of Refs. 49r|52|, the Faddeev equations 53j are solved, having as driving term 



the chiral two-body scattering matrices for the different pairs of the system. In these way, 
the input two-body t matrices in the Faddeev equation contain the information related to 
the generation of the corresponding two-body resonances. 

In line with the above mentioned works, a different strategy to the one discussed in the 
previous section to study the Dj o(980) and D S K systems would be to consider / (980) 
and -D s . (2317) as molecular resonances formed, respectively, in the KK and DK systems 



together with their respective coup 



ed channels and study the three-body system DKK 



following the approach of Refs. 49m2j] . To do this, we consider 10 coupled channels for 



total charge zero and charm C = 1: D°K + K~, D°K°K°, D 
D°TT°r], D + K°K~, D+it-it , D+tt'T}, D + it°it-. 



7T 1 7T 



7T 



D°TT°ir , 



As mentioned above, to solve the Faddeev equations for the DKK system and coupled 
channels, we first need to determine the two-body scattering matrices t for the different 
pairs of the system. This is done by solving the Bethe-Salpeter equation through its on-shell 
factorization form 



20 



24 



26 



36] 



t = (i - vgy l v. 



(28) 



The kernel V in Eq. (1281) corresponds to the lowest order two-body amplitude obtained from 
a suitable Lagrangian and Q represents the loop function of two hadrons. 

In case of the DKK system and coupled channels, we have two different types of inter- 
actions: one involving two light pseudoscalars, like KK, nir, and the other between a heavy 
and a light pseudoscalar meson, like DK, Dtt. 



For the description of the KK system, we follow Refs. 



24 



25] and solve Eq. ([28]) consid- 
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ering KK, tttt and ttt] as coupled channels. The kernel V is obtained from the lowest order 
chiral Lagrangian for the process PP — > PP, with P representing a light pseudoscalar (i.e., 

7T, K, 7]) 



C 



PP 



12f 



■{(d^PP-Pd^Pf + MP 4 ). 



(29) 



In Eq. (l29l) . / is the pion decay constant, the symbol () indicates the trace in the flavor 
space of the SU(3) matrices appearing in P, P is a matrix containing the different Goldstone 
bosons and M a mass matrix: 



7T 



(30) 



M 



/ 



m 2 
m 2 









2m 2 K 



m 



(31) 



The V matrix obtained using the Lagrangian of Eq. (1291) is a function of the Mandelstam 
variables s, t and u. This matrix is further projected on s-wave and the resulting expressions 



can be found in Ref. 



24|. 



The 
of Ref. 



oop function Q in Eq. ( 128]) is calculated using the dimensional regularization scheme 



25j. In the present case, i.e., for a two pseudoscalar system, 



16n 2 



, Or 

E 



a r (fi) + In— f 4 



m lr m 2r ~ m lr + ^ 



2E 2 



In 



m 



2r 



711 



It 



In (£ 2 - (m 2 r . - ml) + 2g r £) + In (e 2 + (m 2 r - m 2 r ) + 2g r £ 



In ( - E 2 + (m 2 r - m\ r ) + 2g r £) - In ( - £ 2 - (m 2 r - m\ r ) + 2g r £ 



(32) 



In Eq. (13"2"|) . is the total energy of the two-body system, m lr , m 2r and q r correspond, 
respectively, to the masses and the center of mass momentum of the two pseudoscalars 
present in the rth channel, /i is a regularization scale and a r (/i) a subtraction constant. 
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Following Ref. 25j, we have taken /x = 1224 MeV and a value for a r (/i) ~ —1 (note that 
there is only one independent parameter here since a change in \i can be reabsorbed in 
a r ). In this way we can reproduce the observed two-body phase shifts and inelasticities for 
the different coupled channels as done in Refs. 24j, |25j. The resulting scattering matrix t 
exhibits poles on the unphysical sheet which are related to the resonances <r(600), /o(980), 
a (980). 

In the case of the subsystem constituted by a heavy and a light pseudoscalar mesons, like 
DK, Dir, since the heavy mesons contain both light and heavy quarks, one expects both the 
chiral symmetry of the light quarks and the symmetry of the heavy quarks to be considered. 
Having this in mind, to determine the scattering matrix of a heavy meson H with a light 



pseudoscalar P, we follow Refs. 



39 



42j, where the leading order Lagrangian describing this 



interaction is given by the kinetic and mass term of the heavy mesons (chiraly coupled to 
pions), 



C = D^HD^H ] - M 2 H HH\ 



(33) 



with H = ( D° D + j collecting the heavy mesons, whose mass in the chiral limit is 
Mhi P is given by Eq. fl30]) and D M is the covariant derivative [48 



= (<9 M + r M )#t, 

D,H = HCd ll + Tl), 



T M = ~{u ] d^u + ud^u ] ) 
u 2 = e^ p l f . 



(34) 



For the process which concerns us, i.e., HP — > HP, Eq. ( 1331) becomes 
C hp = {d»H[P, d,P]H^ - H[P, d,P}d^H^} , 



(35) 



and the lowest order amplitude obtained from this Lagrangian in terms of the Mandelstam 
variables reads as 



Vis 



U 



(36) 



In Eq. (1361) i and j represents the initial and final channels, res pec tively, and the CV,- co 



efficients have been earlier calculated and can be found in Refs. 
further projected on s-wave. 



39 



42[. This potential is 
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As in Ref. 39|, |42|, we consider the coupled channels DK, D s r] and D s n for strangeness 
-1 and Dtt, Drj, D S K for strangeness 0. The loop function Q of Eq. f )28|) is determined 



using Eq. (132 p with /i = 1000 MeV and a = —1.846 [41|, |42[, obtaining in this way a pole in 
the DK system for total isospin at 2318 MeV, which corresponds to the state D s * (2317), 
and a pole at 2446-i43 MeV in the Dtt system in isospin 0, associated with the resonance 
£>*(2400). 

Once the two-body scattering matrices are calculated, we can proceed with the determi- 



nation o 
of Refs. 



' the three-body T matrix for the DKK system. To do that we use the approach 
in which the Faddeev partitions, T 1 , T 2 and T 3 , are written as 

3 



1,2,3, 



(37) 



with ki (k'j) being the initial (final) momentum of the particle i and f the two-body t-matrix 
which describes the interaction of the (jk) pair of the system, j ^ k ^ % — 1, 2, 3. The total 
three-body T-matrix is obtained by summing the T l partitions, 

3 

T 2 + T 3 = ^ 53 (^ " + T R (38) 



T = T 



where we define 



3 3 
i=l j^i=l 



(39) 



The T l jl partitions in Eq. (!37|) satisfy the following set of coupled equations 



Tg = fg'H 3 + f \G l3i T£ + G ijk Tl k ], i ^ j, j ^ k = 1, 2, 3. 



(40) 



where corresponds to the three-body Green's function of the system and its elements are 
defined as 



1 



(41) 



2E k (k'i + kj) JV~s~ Ei(kQ - Ejikj) - E k (k[ + kj) + ie 

with iVjfc = 1 for mesons and Ei, 1 — 1,2, 3, is the energy of the particle /. 

The G l i k matrix in Eq. (140]) represents a loop function of three-particles and it is written 

as 



Qij k 



d 3 k" 

W) 
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(42) 
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with the elements of being 

~gHk»,s lm ) = -2^^ J: , 

2E l {k") 2E m (k") ^ - E x {k") - E m {k") + xe 

and the matrix E l ^ h , with explicit variable dependence, is given by 



(43) 



E^ k (k", k' k k , O = ti{s k Z)<?\k"i h) W\k' h) 



tHsr 



(44) 



In Eq. ( l4"3"j) . Jsi m is the invariant mass of the (Im) pair and can be calculated in terms of 
the external variables. The upper index k" for the invariant mass s^ u of Eq. (j4"4"j) indicates 
its dependence on the loop variable (see Ref. 50| for more details). 

The T]j partitions given in Eq. fj4*0|) are functions of the total three-body energy, y^i, 
and the invariant mass of the particles 2 and 3, y/s^. The other invariant masses, yfsy). 



and y/s^i can be obtained in terms of y/s and ^/s23, as it was shown in Ref. [50|, |5l|. In 
this model, peaks obtained in the modulus squared of the three-body T-matrix are related 
to dynamically generated resonances. Finally, it should be mentioned that the first term 
in Eq. ( 138]) can not give rise to any state generated due to the three-body dynamics and, 
hence, we can just study the properties of the Tr matrix defined in Eq. ( |39|) . 




V^T(MeV) 



V s (MeV) 3040 



FIG. 4: Squared amplitude for the DKK channel for total isospin 1=1/2 with the KK 

subsystem in isospin zero. 



Further, we work in the charge basis, then, to associate the peaks found in the three- 
body T-matrix with physical states we need to project T on an isospin basis. We do this by 
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defining a basis where the states are labeled in terms of the total isospin I of the three-body 
system and the isospin of one of the two-body subsystems, which in the present case is 
taken as the isospin of the KK subsystem or the system made by particles 2 and 3, J23, 
and evaluate the transition amplitude (I, I23\Tr\I, I23). The isospin 7 2 3 can be either or 
1, thus, the total isospin I can be 1/2, or 3/2. For the cases involving the states \I = 1/2, 
I23 = 1), \I = 3/2, J23 = 1) we find no structure which could be related to a resonance or a 
bound state. Thus, in the following, we discuss the case 1 = 1/2 with I23 = 0, where we do 
find a resonance. 

In Fig. H]we show the results obtained for the modulus squared of the scattering amplitude 
of the DKK channel for total isospin 1/2 with the KK system in total isospin 0. A peak 
around 2890 MeV with a width of 55 MeV is found when the KK system is in the isospin 
zero configuration with an invariant mass of around 985 MeV, thus, forming the / (980) 
resonance. 

If instead of using the isospin base \I, I23) we use \I, 112), with ii 2 the isospin of the (12) 
subsystem, i.e., the DK system, no clear signal for the peak shown in Fig. H] is observed 
when the DK system is in isospin zero. 

These results are very similar to the ones obtained in the previous section using QCD 
sum rules, in which a state of mass around ~ 2900 MeV is found to couple strongly to the 
£>/o(980) system. 

IV. SUMMARY 

We have studied the DKK system using two different methods: one based on QCD 
sum rules and other on solving few-body equations. In the former case, the Df (980) and 
D s *o(2317)K configurations of the DKK system have been investigated and a state with 
a mass 2.88±0.36 GeV which couples strongly to a Dfo molecular current is obtained. In 
the latter, the Faddeev equations have been solved with input two-body t matrices which 
generate the /o(980) and _D s *o(2317), respectively, in the KK and DK systems and related 
coupled channels. As a result, a state with a mass close to 2.9 GeV and a width of 55 MeV 
was found when the KK subsystem generates the /o(980) resonance. The findings obtained 
within these two different methods are in fact very similar, hinting towards the existence of a 
£>/o(980) molecular state with a mass close to 2.9 GeV. A state with this mass has not been 

18 



discovered experimentally so far, and the heaviest known D meson is the .0(2750), whose 
mass is around 150 MeV below the one found in this manuscript. We strongly encourage 
the search of a state decaying into DKK with the characteristics of the one found here. 
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